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ABSTRACT 

We study the complex structure of the Bullet cluster radio halo to determine the Dark 
Matter (DM) contribution to the emission observed in the different subhalos corre¬ 
sponding to the DM and baryonic dominated regions. We use different non-thermal 
models to study the different regions, and we compare our results with the available 
observations in the radio, X-ray and gamma-ray bands, and the Sunyaev-Zel’dovich 
(SZ) effect data. We find that the radio emission coming from the main DM subhalo 
can be produced by secondary electrons produced by DM annihilations. In this sce¬ 
nario there are however some open issues, like the difficulty to explain the observed 
flux at 8.8 GHz, the high value of the required annihilation cross section, and the lack 
of observed emission coming from the minor DM subhalo. We also find that part of 
the radio emission originated by DM annihilation could be associated with a slightly 
extended radio source present near the main DM subhalo. Regarding the baryonic 
subhalos, the radio measurements do not allow to discriminate between a primary or 
secondary origin of the electrons, while the SZ effect data point towards a primary 
origin for the non-thermal electrons in the Main Subcluster. We conclude that in or¬ 
der to better constrain the properties of the DM subhalos, it is important to perform 
detailed measurements of the radio emission in the regions where the DM halos have 
their peaks, and that the separation of the complex radio halo in different subhalos is 
a promising technique to understand the properties of each specific subhalo. 
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1 INTRODUCTION 

Several observations suggest that most of the matter in the 
universe is dark and non-baryonic (see, e.g., Bertone et al. 
2005 for a review, and references therein). Various methods 
based on direct or indirect detection have been proposed to 
constrain the nature and the properties of the Dark Mat¬ 
ter (DM). A promising way to study the DM properties 
consists in the detection of the electromagnetic emissions 
coming from astrophysical objects as a consequence of DM 
particles annihilation and secondary particle production, 
because the spectral features of these emissions are expected 
to be closely correlated with the nature, the composition and 
the mass of the DM particles (see, e.g., Colafrancesco et al. 
2006). 

A recent study (Colafrancesco et al. 2015) showed that 
many cosmic structures with different size and mass can be 
promising candidates to detect the DM signal, especially 
in the radio band, using future high-sensitivity experiments 

* E-mail: Paolo.Marchegiani@wits.ac.za 
f E-mail: Sergio.Colafrancesco@wits.ac.za 


like the Square Kilometre Array (SKA) (e.g., Dewdney et 
al. 2012). Among these structures, dwarf spheroidal galaxies 
appear to be good candidates because in these objects the 
DM component is supposed to be dominant with respect to 
other diffuse components of baryonic origin (e.g., Regis et al. 
2014), whereas in galaxy clusters the DM signal is expected 
to be stronger, but other mechanisms of baryonic origin can 
produce emissions superposed to the DM emission (see, e.g., 
Colafrancesco et al. 2010 for the case of the Perseus cluster). 

Non-thermal emission in the atmospheres of galaxy 
clusters manifests itself with a wide range of emission mech¬ 
anisms and can be then observed over a wide range of fre¬ 
quencies. Synchrotron emission from electrons either primar¬ 
ily accelerated at shocks (Tribble 1993) or re-accelerated by 
turbulence (Brunetti et al. 2001), or secondarily produced 
by proton-proton collisions (Blasi & Colafrancesco 1999) or 
DM annihilation (Colafrancesco & Mele 2001, Colafrancesco 
et al. 2006, 2011b) is expected to be observed in the ra¬ 
dio band. The same primary and/or secondary electrons 
are also expected to emit by inverse Compton Scattering 
(ICS) off the Cosmic Microwave Background (CMB) pho¬ 
tons (and possibly other cosmic backgrounds), and by non- 
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thermal bremsstrahlung off the nuclei of the thermal gas 
over a wide energy range from nltraviolet to gamma-rays 
(e.g. Blasi & Colafrancesco 1999). A strong gamma ray emis¬ 
sion, due to the P —>■ 77 decay, is expected in secondary 
models, either produced by p — p collisions or by DM anni¬ 
hilation (e.g., Colafrancesco & Blasi 1998, Colafrancesco et 
al. 2006). Other astrophysical effects related to these non- 
thermal components are also expected, like a non-thermal 
Sunyaev-Zel’dovich effect (SZE) (e.g. Colafrancesco et al. 
2003) and gas heating by Coulomb and hadronic interac¬ 
tions (e.g. Colafrancesco & Marchegiani 2008). 

For these reasons, the detection of genuine DM- 
produced signal in galaxy clusters appears to be not easy, 
especially when the baryonic and the DM emissions are co- 
spatially located. Therefore, a good way to search for the 
DM emission is to study galaxy clusters where the DM and 
the baryonic halos are spatially offset. The most famous clus¬ 
ter with these properties is the Bullet cluster (IE 0657-558), 
where observations of the X-ray emission (e.g. Markevitch 
et al. 2002) and of the gravitational leasing (Clowe et al. 
2006, Bradac et al. 2006) showed that distributions of the 
hot gas in the Intra Cluster Medium (ICM) and of the DM 
are centered at different positions. The other case of A520 
where there is some spatial offset between DM and baryons 
is not a very clean case to study (Mahdavi et al. 2007, Clowe 
et al. 2012 ). 

The Bullet cluster is a merging cluster which shows a 
complex structure and, because of the spatial offset between 
hot baryons, DM and galaxy distributions, appears to be 
one of the best laboratories to decipher the stratification of 
the complex electron population. In the X-rays this cluster 
shows the presence of two subclusters: a Main Subcluster 
(MS), with an elliptical shape elongated in the North-South 
direction, and a Bullet Subcluster (BS), located on the west 
of the MS, with a more asymmetric shape that shows a front 
shock in the west direction. This is currently interpreted by 
assuming that the cluster is in a post-merging phase, and 
the BS has just crossed the MS, producing gas heating and 
particles acceleration. The gravitational leasing analysis of 
the Bullet cluster shows two main halos: the larger one is lo¬ 
cated on the east of the MS, and the smaller one is located 
on the west of the BS. The offset between the hot gas and 
DM halos is explained assuming that this is a consequence 
of the recent merging where the collisional (the baryons) 
and the collisionless (DM) components have been separated 
after the powerful merging event. Finally, the spatial dis¬ 
tribution of the galaxies in the cluster seems to be related 
more with the DM distribution rather than with the hot gas 
distribution. 

Radio observations of the Bullet cluster (Liang et al. 
2000, Shimwell et al. 2014) show a very complex morphology, 
with an elliptical radio halo elongated in the east-west direc¬ 
tion, a large radio relic on the east of the halo (Shimwell et 
al. 2015), and several mini halos around galaxies and point¬ 
like sources. Interestingly, the overall radio halo, after the 
sources removal procedure (see Shimwell et al. 2014), ap¬ 
pears to have three different peaks, two of which are located 
near the X-rays peaks, and another one located near the DM 
Eastern (DME) region. The spectrum of the radio halo also 
shows very interesting features: while the emission from the 
internal part of the halo, where the MS is located, has an al¬ 
most perfect power-law shape, the spectrum of the extended 


region shows a flattening of the spectrum at ly ^ 2.3 GHz, 
and a steepening at rz > 5 GHz. This fact can be produced 
by the presence of different emission components with differ¬ 
ent spectra that dominate alternately in different frequency 
ranges, one of which can be due to the DM component. 

For this reason, it is interesting to study in detail the 
emissions coming from the different regions of the radio halo 
and explore the possibility to separate the baryonic and the 
DM-produced emissions. This goal requires to proceed with 
a very detailed study of the different halos, using appropriate 
models for each of them, starting from radio emission and 
using, when available, information obtained in other spec¬ 
tral bands. 

In fact, due to the degeneracy between electrons density 
and magnetic field intensity arising from radio observations 
(e.g., Longair 1994, Colafrancesco et al. 2005), it is impor¬ 
tant to observe, or constrain, the non-thermal emission in 
the Bullet cluster in other spectral bands, like the X-ray and 
gamma-ray bands. Unfortunately, only upper limits have 
been obtained so far. The most recent NuSTAR observa¬ 
tions in the hard X-rays band obtained with an integration 
time of 266 ks (Wik et al. 2014) provided an upper limit 
E(50 — 100 keV) < 1.1 x 10“^^ erg cm“^ s“^. In the gamma 
ray energy range, an upper limit has been derived by the 
Fermi-LAT experiment, 7^(0.2 - 100 GeV) < 2.75 x 10"® 
cm"® s"^ (Ackermann et al. 2010). The same authors pro¬ 
vided in their Figure 1 a gamma ray upper limit, extrapo¬ 
lated to E > 0.1 GeV, that is E(> 0.1 GeV) < 5.7 x 10"® 
cm"® s"^. 

Another source of information on the properties of the 
Bullet Cluster is provided by the SZE observations; the 
study of its spectral shape, especially at high frequencies 
(ly > 300 GHz), where the relativistic corrections are more 
relevant (Colafrancesco et al. 2003, Colafrancesco et al. 
2011 a), allows to obtain information on the thermal and pos¬ 
sibly the non-thermal electrons in the cluster atmosphere. 
In the Bullet cluster, the data at high frequencies obtained 
with Herschel-SPIRE (Zemcov et al. 2010) allowed to es¬ 
tablish that, in addition to main thermal component, there 
is a second electron population (more likely of non-thermal 
origin) that is contributing to the total observed SZE (Co¬ 
lafrancesco et al. 2011 a), and that the gas temperature dis¬ 
tribution along the line of sight is strongly inhomogeneous 
(Prokhorov & Colafrancesco 2012). 

By using all this information, we perform in this paper 
a detailed study of the Bullet cluster, starting from the spec¬ 
tral and the spatial properties of the radio halo, and using 
different emission models that are appropriate for the dif¬ 
ferent regions. We also check the information obtained from 
the radio band considering the corresponding emissions in 
the high-energy (i.e.. X-ray and gamma ray) bands, com¬ 
paring the predictions with the available upper limits and 
with the sensitivities of coming and planned instruments. 
We also study the SZE in a more detailed way than in pre¬ 
vious studies, and make predictions for future observations 
with instruments like Millimetron, a planned cooled 10 m 
space telescope which will operate in the Lagrangian L2 
point in the 20 pm - 20 mm wavelength range (Smirnov 
et al. 2012; Kardashev et al. 2015). 

The outline of the paper is the following: in Sect. 2 we 
study the structure of the radio halo emission in the Bul¬ 
let cluster with specific interest to the various sub-peaks of 
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this emission. We discuss in Sect.3 the consequences of the 
radio halo analysis on the DM and CR distribution models, 
the possible issues, and the predictions in the high-energy 
bands. We then focus in Sect.4 on the analysis of the SZE. 
We finally summarize our results and draw our conclusions 
in Sect.5. 

Throughout the paper, we use a fiat, vacuum-dominated 
cosmological model following the results of Planck, with 
Dm = 0.308, Da = 0.692 and Hq = 67.8 km s“^ Mpc“^ 
(Ade et al. 2015). With these values the luminosity distance 
for the Bullet cluster at 2 ; = 0.296 is Dl = 1576 Mpc, and 
1 arcmin corresponds to 273 kpc at this distance. 


2 THE STRUCTURE OF THE RADIO HALO 

The observed radio halo of the Bullet cluster has an ellip¬ 
tical shape with the major axis oriented in the East-West 
direction (see, e.g., Fig.5 in Shimwell et al. 2014). Actually, 
the radio brightness inner contours of the radio halo suggest 
that its overall elliptical shape can be due to the combi¬ 
nation of emissions of different subhalos; in particular, the 
main peak is located near the center of the MS X-ray emis¬ 
sion, and the radio emission appears to be elongated both in 
the west direction, with a secondary peak located near the 
BS X-ray peak, and in the south-east direction, with another 
peak located near the peak of the DME distribution. A little 
residual emission seems to be present also near some of the 
radio sources identified in the Bullet cluster, as the sources 
indicated with A, G, H, F and K in Shimwell et al. (2014). 

It is therefore justified to test the hypothesis that the 
overall radio emission is mainly due to the combination of 
three distinct radio sub-halos. In this framework we will first 
test some hypotheses on the origin of the radio halo emis¬ 
sion that are consistent with both the available observations 
and the theoretical expectations, and then we will verify the 
consequences of these hypotheses based on their consistency 
with the available observations. 

The hypothesis on the origin of the Bullet cluster radio 
halo that we want to test is the following: 

i) the sub-halos coincident with the MS and the BS are of 
baryonic origin; we assume that they are produced by a sec¬ 
ondary electron model (SEM) of hadronic origin where the 
CR protons are normalized according to the Warming Ray 
(WR) model (see Colafrancesco & Marchegiani 2008 for de¬ 
tails) that is able to reproduce the X-ray temperature distri¬ 
bution of the two baryonic halos. The use of SEM seems to 
be reasonable because in this model the strongest emission 
comes from the regions where the ICM density is highest; 
the normalization of the non-thermal protons following the 
WR model provides actually an upper limit on their density, 
because a larger density would provide a stronger gas heat¬ 
ing and therefore a higher temperature than the observed 
one; 

ii) the radio sub-halos coincident with the DM halos are pro¬ 
duced by secondary electrons produced by DM annihilation; 
we will focus on the case where the DM is constituted by 
neutralinos, i.e. the lightest particles in the minimal super- 
symmetric extension of the standard model (MSSM). 


2.1 The Main Subcluster 


The main subcluster (MS) is found at the location where the 
X-ray emission has its main peak (Markevitch et al. 2002), 
and at the same location also the radio brightness maps 
show a peak coincident with the X-ray one (Shimwell et al. 
2014). 

The thermal gas temperature in this region is kT = 
14.2lto’2 keV (Wik et al. 2014), and the gas density profile 
can be approximated by an isothermal / 1 -model: 


nth{r) = no 


1 + 



( 1 ) 


with no = 7.1 x 10~® cm~®, /3 = 1.04 and rc = 516 kpc 
(Ota & Mitsuda 2004), where the approximation of spheri¬ 
cal symmetry is used. The spatial extent of the radio emis¬ 
sion in this region corresponds to the smallest of the regions 
considered in the analysis of Liang et al. (2000) (region ^ 
3 in their Figure 4), and the relative radio spectrum is the 
lower flux spectrum plotted in their Figure 7, which shows 
a radio spectral index of an ~ 1.4. 

When we apply the WR model to this subcluster, in 
which we assume that the spatial profile of non-thermal pro¬ 
tons is oc nth{r)°‘, we find that this model requires a = 1 , 
consistently with many other isothermal clusters (see the 
analysis of Colafrancesco & Marchegiani 2008 for details). 
Assuming a proton spectrum of the type 


Nph) = 


( 2 ) 


with Sp = 2.7, the WR model provides the value Np^o = 4.1 x 
10“® cm“® for the CR proton normalization, which implies 
a ratio between the CRs and the thermal gas pressures at 
the cluster center of Pc r/P th = 0.42. 

The radius of the X-ray emission in this region is 
R ~ 520 kpc, but, to determine the value of the magnetic 
field necessary to produce the observed radio emission, we 
calculate the emission produced inside a radius of ~ 160 kpc, 
from the size of region # 3 of Liang et al. (2000). In this way, 
if we assume that the magnetic field is uniform through¬ 
out this region, we find that the diffuse radio emission is 
reproduced by secondary electrons emitting in a magnetic 
field with amplitude B = 7.5 /rC. Such value is compati¬ 
ble with the overall findings in other clusters (e.g., Carilli & 
Taylor 2002), and it is reasonable in a strong merging clus¬ 
ter like the Bullet cluster which is very hot, with a highly 
non-uniform gas distribution where the recent shock activ¬ 
ity could have amplified the pre-existing magnetic field (see 
also Shimwell et al. 2015). Figure [T] shows that this model 
can indeed reproduce quite well the data observed in the 
central region studied by Liang et al. (2000). 

The hard X-ray flux produced in the 50-100 keV band 
by the secondary electrons in the WR model is E(50 — 
100 keV) = 1.8 X 10“^® erg cm“^ s“^, much lower than 
the NuSTAR limit (Wik et al. 2014). The gamma ray flux 
at E > 0.1 GeV, mainly produced by the vr® —>■ 77 channel, 
expected in this model is F{E >0.1 GeV) = 5.1 x 10“^® 
cm“^ s“^, well below the Fermi-LAT limit reported by Ak- 
ermann et al. ( 2010 ). 

This model for the MS is therefore fully consistent with 
all multi-frequency data available for the cluster, from radio 
to gamma-rays. 
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Figure 1. The radio emission in the region of the Main Subcluster 
fitted with a WR model with Sp = 2.7 and B = 7.5 fiG. The data 
shown in this plot correspond to the central region of the cluster 
as given by Liang et al. (2000). 

2.2 The Bullet Subcluster 

The BS is centered on the region of the second peak in the 
X-ray emission, i.e. the one that corresponds to the “bullet” 
that crossed the MS and is producing the shock front. 

This subcluster is colder in its central part (with a tem¬ 
perature of kT ~ 7 keV), and it is hotter at its periphery 
(where the temperature reaches kT ~ 18 keV), according to 
the analysis of Markevitch et al. (2002). The same authors 
found that the gas density is nth ~ 1-7 x 10“^ cm“® in the 
cold region, within a radius of ~ 50 kpc from the center of 
the BS, and is nth ~ 4.1 x 10”^ cm“^ in the hotter region, 
within a radius of ~ 260 kpc from the BS center; in both re¬ 
gions it is found that the density is approximately constant 
with radius. We can therefore consider that the BS has a 
structure similar to a cool-core cluster. 

When we apply the WR model to the BS we obtain a 
value a = 0.39. For a value of the protons spectral index 
Sp — 2.7, the normalization of the proton spectrum required 
by the model is Np^ = 4.9 x 10“® cm“® {PcR/Pth = 0.46), 
while for Sp = 2.9 we find Npfi = 5.8 x 10“® cm“® 
{PcR/Pth = 0.67). 

We use this model to describe the radio halo data in the 
region # 2 in Figure 4 of Liang et al. (2000), which is similar 
in extension to the smaller region considered in the analy¬ 
sis of Shimwell et al. (2014), and hence we use the higher 
level radio flux data in the spectrum obtained by Liang et 
al. (2000), i.e. those in their Figure 7, and the lower ampli¬ 
tude radio flux spectrum reported in Figure 7 of Shimwell et 
al. (2014). The region ^ 2 in Liang et al. (2000) is actually 
larger than the BS, and includes also the DME region and 
some radio sources. So, now we compare the radio emission 
produced in the BS with the region ^ 2 data, but in the fol¬ 
lowing we will consider the contributions from other sources 
in this region. The data sets of Liang et al. (2000) and those 
of Shimwell et al. (2014) show that the radio halo spectrum 
in this region has a power-law shape with index an ~ 1.6 
up to a frequency of v ~ 2.3 GHz, while at higher frequen¬ 
cies the spectrum flattens between 2.3 and 5 GHz with a 
further steepening at frequencies v ^ 5 GHz. This spectral 
shape might suggest that the observed spectrum is due to 


Figure 2. The radio emission in the region of the Bullet Sub¬ 
cluster fitted with two WR models with Sp = 2.9 and B = 60 fiG 
(solid line) and with Sp = 2.7 and B = 27 fiG (dashed line). The 
data shown in this plot correspond to the region # 2 of the clus¬ 
ter as given by Liang et al. (2000) (triangles) and to the smaller 
region considered by Shimwell et al. (2014) (asterisks). 

the superposition of different components, one of which can 
be the BS, and a second one is flatter and dominates the 
overall spectrum at 2.3 m 5 GHz, and then steepens at 
higher frequencies. 

In Fig[2]we show the radio spectrum produced in the BS 
for the two considered WR models, with a magnetic held of 
B = 60 /iG for Sp = 2.9, and with B = 27 /rG for Sp = 2.7. 
These values of the magnetic field are clearly extreme, so 
probably the BS emission alone is not sufficient to produce 
the flux observed in the whole region # 2 of Liang et al. 
(2000); anyway, the BS is located near the shock region, 
and it hosts a cool core, hence a relatively high value of 
the magnetic field is conceivable. As we can see, the first 
model reproduces well the data up to a frequency of 2.3 GHz, 
and fails to reproduce the data at higher frequencies, while 
the second model is less accurate at small frequencies but, 
given the large data dispersion, it can explain the observed 
spectrum at all the frequencies. 

The IGS flux produced in the 50-100 keV band is 4.1 x 
10“^® erg cm“^ s“^ for Sp = 2.9, and 3.9 x 10“^® erg cm“^ 
s“^ for Sp = 2.7, both much lower than the limit of Wik et 
al. (2014). 

The gamma ray flux for B > 0.1 GeV produced by P 
decay for Sp = 2.7 is 9.7 x 10“^^ cm“^ s“^, and for Sp = 2.9 
is 9.6 x lO”^*^ cm“^ s”*^, both below the limit of Ackermann 
et al. (2010). 


2.3 The Dark Matter Eastern peak 

Gravitational leasing image of the gravitational potential 
of the Bullet cluster shows two peaks separated from the 
X-ray peaks (Clowe et al. 2006), that are supposed to be 
dominated by the Dark Matter density distribution. The 
luminous component (i.e., the galaxies) seems to follow the 
DM distribution rather than the thermal gas distribution. 
The DME peak is found to be the most massive one, with 
a mass Mdme ~ 5.6 x within a region of 520 kpc 

radius (see Bradac et al. 2006). 
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Figure 3. Annihilation time (as defined in eq|4jl for a DM halo 
with M-^ = 500 GeV and the density parameters of the DME 
region for (av) = 4.3 X 10“^^ cm® s“^ (solid line), (av) = 3 X 
10“^® cm® (dashed line) and (av) = 3 X cm® s”^ 

(dot-dashed line), compared with the Hubble time at the Bullet 
Cluster redshift (long dashed line). 


Following the procedure described in Colafrancesco et 
al. (2015), we study the case where the DM is constituted 
by neutralinos with different values of the mass and compo¬ 
sitions, and we use the relations between the halo mass and 
the parameters of the DM distribution, as the concentration 
parameter and the central density, and, for a Navarro, Frenk 
& White (1996) mass density distribution. 


p{r) 




(3) 


we obtain the values Ps =6.14 x 10®Pc and = 378.5 kpc, 
where pc is the critical density of the universe. 

Formally, this DM density profile tends to infinite value 
at r —>■ 0; indeed, this is not physically possible because the 
increasing DM density at the cluster center increases the 
annihilation probability, and therefore the annihilation time 
scale 


iann R ) 


1 

nDM{r){av) ’ 


(4) 


where udmIt) = p{r)/M^ is the neutralino numerical den¬ 
sity, and (av) is the mean value of the annihilation cross sec¬ 
tion, decreases. When the annihilation time given in eq. O 
is lower than the Hubble time tn at the cluster redshift, 
the annihilation processes are dominating and, as a conse¬ 
quence, the DM density cannot increase again. So, one can 
find a minimum radius over which the DM distribution has 
the shape of eq. m, and under which the density is roughly 
constant. In Figure [3] we show the annihilation time for a 
neutralino mass of = 500 GeV and for three values of 
the cross section: a reference value of (av) = 3 x 10“^® cm® 
s“^ (this is an estimate of the DM relic abundance value), 
a value of the order of the cross section we will find in next 
analysis, (av) = 4.3 x 10“^^ cm® s“^, and an intermediate 
value of (av) = 3 x 10“^^ cm® s“^. As we can see, the ra¬ 
dius at which the DM density stops raising is quite small; 
for a cross section value of (av) = 4.3 x 10~^^ cm® s“^ it is 
~ 3 x 10“® pc, and it is smaller for decreasing cross section 
values. 


By using this radial distribution density for DM, we 
calculate the production rate of secondary electrons and the 
equilibrium spectrum due to energy losses following the pro¬ 
cedure described in Colafrancesco et al. (2006). In this cal¬ 
culation, we assume different masses and compositions for 
the neutralino by using the DarkSusy package (Gondolo et 
al. 2004). We assume here a boosting factor due to the sub¬ 
halos DM distribution (Fieri et al. 2011) equal to 1, and 
larger values can be considered eventually simply by multi¬ 
plying the result by this factor (see, e.g., Colafrancesco et 
al. 2011b). 

First, we check if the spectral shape of the radio emis¬ 
sion from DME can explain the flattening of the radio spec¬ 
trum at 1 / (5 2.3 GHz. In Figure 0] we show the emissions 
produced using the WR model in the MS and BS regions, 
using for the BS the model with Sp = 2.9, that reproduces 
well the data up to 2.3 GHz, and we show the radio spectrum 
produced in the DME region for = 500 GeV, composi¬ 
tion W'^W~, B = 0.015 pG and (av) = 8.1 x 10“^^ cm® 
s“^. We discuss this model because it produces the best re¬ 
sult, regarding to the spectral shape of radio flux w.r.t. the 
observed one, among all the models we considered, for neu¬ 
tralino masses of — 9, 60 and 500 GeV, compositions 
bb, t'^t~ and W'^W~, and by searching for the value of the 
magnetic field that produces the spectral shape closest to 
the observed one, and calculating the corresponding cross 
section to fit the radio data. The radio emission from DME 
is subdominant for v ^ 2.3 GHz, and is dominant at higher 
frequencies reproducing the correct slope up to 5.9 GHz, 
whereas it overestimates the point at 8.8 GHz. In addition 
to this problem, this model is clearly unrealistic, because it 
requires a value of the magnetic field smaller than the lower 
limit of 0.2 p,G found with HXR measurements by Wik et 
al. (2014), and a DM annihilation cross section larger by 
a factor of ~ 10® w.r.t. the upper limit (av) ~ 10“^“^ cm® 
s“^ obtained by the Fermi analysis of dwarf galaxies (Acker- 
mann et al. 2014). Anyway, it is possible to obtain a similar 
spectrum by changing the magnetic field value and the cor¬ 
responding cross section (see Figure [5)|. For higher values of 
the magnetic field, the resulting spectrum shows a smaller 
curvature at high frequencies, so it overestimates the point 
at 8.8 GHz by a larger factor w.r.t. the previous model, but 
it is compatible with the HXR measures for B ^ 0.055 pG 
and, for 5^5 pG, it requires a cross section higher than 
the Fermi upper limits by a factor less than 10®, that can 
be recovered by assuming the presence of a boosting fac¬ 
tor originated by the substructures in the halo. We stress, 
however, that the problem of the spectral steepening at 8.8 
GHz is present in each of these models; we will discuss this 
problem in Sect. 3.1 below. 

For the model with = 500 GeV and composition 
W'^W~, the gamma ray flux produced in the case with B = 
1 pG is 1.8 X 10~®® cm“^ s“®, and for B = 0.015 pG it is 

9.3 X 10“® cm“^ s“® (greater than the Fermi upper limit in 
the last case). 

For comparison, we show also the results for two models 
with = 9 GeV (see Figure O: they reproduce well the 
observed steepening GHz, and provide ICS fluxes of 

4.4 X 10“®® erg cm“^ s“® for the bb model, and 8.0 x 10“®®" 
erg cm“^ s“® for the r^r“ model, both smaller than the 
observed upper limit of Wik et al. (2014). However, these 
models overestimate the radio data at ^ 3 GHz. 
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Figure 4. Radio emission in the different regions of the Bullet 
cluster with: a WR model with Sp = 2.7 and B = 7.5 fiG in 
the MS (solid line), a WR model with Sp = 2.9 and B = 60 fiG 
in the BS (dashed line) and a DM model with = 500 GeV, 
composition W'^W~, B = 0.015 /rG and (av) = 8.1 X 10“^^ cm® 
s“^ in the DME region. Data are taken from region # 3 of Liang 
et al. (2000) (diamonds), from region # 2 of Liang et al. (2000) 
(triangles) and from the smaller region of Shimwell et al. (2014) 
(asterisks). 
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Figure 5. Radio emission in the different regions of the Bullet 
cluster with: a WR model with Sp = 2.9 and B = 60 fiG in the 
BS (solid line) and several DM models with = 500 GeV and 
composition W^W~ in the DME region with B = 0.015 fiG and 
(av) = 8.1 X 10“®^ cm® s“® (long dashed line), B = 1 fiG and 
{av) = 1.6 X 10“^'^ cm® s“® (three dotted-dashed line), B = 5 
fiG and (av) = 1.1 X 10“^^ cm® s“^ (dot-dashed line), B = 10 
fiG and (av) = 6.5 X 10“^^ cm® s“® (dashed line), B = 20 /rG 
and (av) = 5.8 X 10“^^ cm® s“® (dotted line). Data are taken 
from region # 2 of Liang et al. (2000) (triangles) and from the 
smaller region of Shimwell et al. (2014) (asterisks). 


2.4 The Dark Matter Western peak 

Now we calculate the radio emission produced by DM mod¬ 
els in the DMW region, where a strong radio component 
does not seem to be present. 

From gravitational lensing results, the DMW halo is 
less massive than the DME, with a mass of Mumw ~ 
3.0 X 10^^M 0 inside 310 kpc (Bradac et al. 2006). The NFW 
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Figure 6. Radio emission in the different regions of the Bullet 
cluster with: a WR model with Sp = 2.9 and B = 60 fiG in the BS 
(solid line) and several DM models with = 9 GeV in the DME 
region, with composition bb, B = 30 fiG and (av) = 4.6 X 10“^® 
cm® s“^ (dashed line), and composition t+t“, B = 10 ^iG and 
(av) = 2.6 X 10~^® cm® s“® (dot-dashed line). Data are taken 
from region # 2 of Liang et al. (2000) (triangles) and from the 
smaller region of Shimwell et al. (2014) (asterisks). 

parameters we derive are ps = 7.41 x 10®pc and Vs = 283.5 
kpc. 

For the model with My. = 500 GeV and composition 
W'^W~ we find that, for the same values of the magnetic 
field and the cross section used for the DME, the radio emis¬ 
sion produced in the DMW is smaller by a factor of ~ 0.6 
than the emission in DME. So, this emission is relatively 
low, but in principle not negligible, considering also that in 
this region, near to the shock front, strong magnetic fields 
and recently accelerated relativistic particles can be present. 
We will discuss this problem in Section 3.2 below. 


3 REFINING AND TESTING THE BASIC 
MODEL 

3.1 High frequencies spectrum 

We discuss here the problem related to the excess flux at 8.8 
GHz, that is difficult to be explained in all the models we 
considered. We consider different possible explanations: 
i) our model is simplified because it takes into account only 
the radio emissions coming from the three described regions, 
while the region # 2 in Liang et al. (2000) contains the con¬ 
tributions from a larger area. Other contributions can flatten 
the total emission between 2.3 and 5 GHz, and steepen it at 
8.8 GHz; 

it) the experimental technique to obtain the data is affected 
by some issue: for example, the points of Liang et al. (2000) 
are corrected for the thermal SZE expected at the same fre¬ 
quencies, and these corrections are calculated under the as¬ 
sumption of an isothermal cluster with the density given by 
two beta models centered on the two X-rays peaks, following 
Andreani et al. (1999). However, the assumption of isother- 
mality can be too inaccurate in this very complex cluster 
(e.g., Markevitch et al. 2002, Prokhorov & Golafrancesco 
2012), and also other SZE components can be present, and 
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they can affect the total flux in different ways at different 
frequencies (see Colafrancesco et al. 2011a); 

Hi) an incorrect source subtraction, which is a delicate oper¬ 
ation mostly for extended sources, can alter the total radio 
halo spectrum. Both Liang et al. (2000) and Shimwell et al. 
(2014) report a detailed list of sources in the field. Among 
the listed sources contained inside the region ^ 2 of Liang 
et al. (2000), the most interesting one appears to be the 
source labeled with an A by Liang et al. (2000) and with 
an L by Shimwell et al. (2014), at coordinates (J2000) RA 
06:58:37.9 and Dec -55:57:25. This is the only source inside 
this field with a quite high flux (~ 20 mjy at 1.3 GHz), and 
it is slightly extended compared with the ATCA resolution 
of 2.7” in Shimwell et al. (2014), so the procedure of the flux 
removal can induce some errors in the total flux estimate of 
the halo. 

The A(L) radio source is located near the galaxy 
2MASX J06583806-5557256, in a region rich of galaxies, at 
the coordinates J2000 RA 06:58:38.0 and Dec -55:57:26. It is 
detected in the infrared bands by WISE and in the 2MASS 
survey with Kg, H and J magnitudes 13.4, 14.0 and 15.0, 
respectively, and in the optical band with B magnitude 20.3 
and redshift 2 = 0.297 (Barrena et al. 2002). Another galaxy 
that is very close the previous one is located at coordinates 
J2000 RA 06:58:38.0 and Dec -55:57:23, with B magnitude 
21.0 and redshift 2 = 0.292 (Barrena et al. 2002). The po¬ 
sition of the A(L) source, located close (and likely corre¬ 
sponding) to the SUMSS source at coordinates J2000 RA 
06:58:37.6 and Dec -55:57:24, is between these two galaxies, 
and it is not clear if it is associated with one of them, or it 
is related to a possible interaction. 

We also note that the position of the radio source A(L) 
is very close to the center of the DME region. This fact sug¬ 
gests also the possibility that the radio source A(L) can be 
produced by the DME halo, in the same way we explored in 
the previous Section. Since for a NEW DM density profile 
the radio emission has a sharp central peak and an extended 
halo component, it is possible that this source is the central 
emission of the DME region that has been subtracted (that 
at the distance of the Bullet cluster has a radius of few arc- 
sec). In this case, to reproduce the total emission from DME, 
the DM model should be normalized to the A(L) source spec¬ 
trum, and its emission should be added to the total halo 
emission. 

The spectrum of the A(L) source can be reproduced 
with a DM model with = 500 GeV, W^IV~ composi- 
tion, B = 10 fiG and (av) = 4.3 x 10“^^ cm® s“^. The same 
flux level can be obtained also with a cross section of the or¬ 
der of the maximum value allowed by the Fermi upper limit 
from Dwarf Galaxies ((uu) ~ 10“^"^ cm® s“^, Ackermann et 
al. 2014) and a boosting factor smaller than 10®. We found 
that models with other values of the neutralino mass or dif¬ 
ferent compositions fail to reproduce the correct slope of the 
radio flux, but we note that there is a degeneracy between 
the value of the magnetic field and that of the cross section 
(Colafrancesco et al. 2006, Colafrancesco et al. 2015). 

In Fig[7]we plot the radio spectrum from the region # 
2, reproduced with the two WR models with different Sp val¬ 
ues we discussed before in the BS, compared with the data 
from this region (triangles and asterisks), and the emission 
of the radio source A(L) explained with the DM model with 
the corresponding data (crosses). In the same figure we also 
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Figure 7. Upper panel: radio emission from different regions in 
the Bullet cluster with a WR model with Sp = 2.9 and B = 60 
/iG in the BS (solid line) and a DM model with = 500 GeV, 
composition W'^W~, B = 10 fiG and {av) = 4.3 X 10“^^ cm® 
s“^ (dashed line) in the DME region, and their sum (long dashed 
line). Data are from region # 2 oi Liang et al. (2000) (triangles), 
from the smaller region of Shimwell et al. (2014) (asterisks), and 
from radio source A(L) from Liang et al. (2000) (crosses). With 
the squares we plot the sum of the cluster data and the radio 
source date at the same frequencies (derived with interpolation 
at the frequencies where they are not available in Liang et al. 
2000). Lower panel: as in the upper panel, but for a WR model 
in the BS with Sp = 2.7 and B = 27 iiG. 


show the sum of the models in the two regions compared 
with the sum of the data of the cluster and the radio source 
(squares); we stress that in order to obtain these data at the 
frequencies not included in the Table 2 of Liang et al. (2000) 
we interpolated the data with a power-law model. As we can 
see, the best case seems to take place for the WR model in 
the BS with Sp = 2.7, even if it underestimates the data at 
5 and 6 GHz. The corresponding magnetic field value is 27 
pG; this is a quite high value, but we notice that the BS is 
located near the shock front, and that usually the magnetic 
fields in cool core clusters are stronger than those in other 
clusters (e.g., Carilli & Taylor 2002). 

The model with = 500 GeV, W'^W~ composition, 
magnetic field 10 nG and cross section (av) = 4.3 x 10“^® 
cm® s“® produces an ICS flux in the 50-100 keV band of 
1.9 X 10“^^ erg cm“^ s“®, smaller than the upper limit de¬ 
rived by Wik et al. (2014). 
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The gamma ray flux for E > 0.1 GeV predicted by 
the same model is 6.3 x cm~^ s“^ (given by the sum 

of 7r° decay and ICS and non-thermal bremsstrahlung from 
electrons), smaller than the limit derived by Ackermann et 
al. (2010). 

3.2 Emission from the DMW region 

By considering the model we discussed previously with i) 
a WR model with Sp — 2.7 and B = 27 /rG in the BS to 
explain the spectrum in the region # 2, ii) a DM model 
with = 500 GeV, composition W'^W~, B — 10 fiG and 
{av) = 4.3 X 10“^^ cm® s“^ in the DME region to explain 
the spectrum of the A(L) radio source, and iii) a WR model 
with Sp = 2.7 and B = 7.5 fiG in the MS to explain the 
emission in the internal region, we consider the contribution 
of a DM model coming from the DMW region, by taking 
in account the smaller mass of this halo and with the same 
properties of the neutralino (mass, composition and annihi¬ 
lation cross section) used in the DME model. The spectrum 
resulting from this model is plotted in the Fig[8] The emis¬ 
sion from the DMW region is lower than the emission from 
the DME, but it is not totally negligible. This might be a 
problem because the radio map does not show a strong emis¬ 
sion in this region; there are some radio sources located in 
this region (the A, D, and E sources in Shimwell et al. 2014; 
in particular the A source is located very near to the DMW 
peak), but their flux is small. Also, in this plot we assumed 
that the magnetic field is the same in the DME and DMW 
regions, while in the DMW we can expect that it is higher 
than in the DME because the DMW region is located near 
the shock. However, it is possible that the magnetic field can 
be amplified by the shock compression in the region to the 
west of the shock front, while it can be reduced in the region 
to the east of the shock front, where the DMW is located, 
because of the post-shock dilution. 

Therefore, it seems necessary to perform a deeper study of 
the radio emission from this region to better understand the 
properties of the DMW halo. 

To compare the combined effect of the different compo¬ 
nents with the observed radio map, we calculate the surface 
brightness profiles at 1.5 GHz for the different models shown 
in the FigO We plot these profiles in FigO centering all the 
profiles on the same point. As we can see, the surface bright¬ 
ness produced in the DMW region inside 0.1 arcmin from 
the center is not negligible, and at the distance of 1 ar¬ 
cmin it is of the order of 0.5 mjy arcmin“^, so it should be 
detected with SKA if integrated in a region of a few arcsec 
size, provided that the sensitivity of SKA at this frequency 
is of the order of ~ ^Jy (Dewdney et al. 2012). 

We show also a simulated radio map of the surface 
brightness profiles plotted in the Fig[9] centering the pro¬ 
files on the corresponding positions on the map. The re¬ 
sult is shown in the first panel of the Fig |16l This map can 
be compared to the observed radio map reported in Fig.5 
of Shimwell et al. (2014). In the simulated map the DM 
emissions have narrow and intense peaks, but also extended 
residual emissions; so, it is possible that the peaks have been 
interpreted as point-like (or slightly extended) sources, while 
the residual halo is included in the overall cluster emission. 
This can be the case for the DME, while, as discussed before, 
it seems to be a more difficult expectation for the DMW. We 
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Figure 8. Radio emission from different regions in the Bullet 
cluster with a WR model with Sp = 2.7 and B = 27 fiG in the BS 
(solid line), two DM models with = 500 GeV, composition 
W^W~, B = 10 IJ,G and (crv) = 4.3 X 10“^^ cm® in the 
DME region (dashed line) and in the DMW region (dot-dashed 
line), and a WR model with Sp = 2.7 and B = 7.5 fiG in the MS 
(three dots-dashed line). Data are from region # 2 of Liang et 
al. (2000) (triangles), from the smaller region of Shimwell et al. 
(2014) (asterisks), from radio source A(L) from Liang et al. (2000) 
(crosses), and region ^ 3 of Liang et al. (2000) (diamonds). 
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Figure 9. Surface brightness profiles at 1.5 GHz for: a WR model 
with Sp = 2.7 and B = 7.5 /rG in the MS (solid line), a WR model 
with Sp = 2.7 and B = 27 iiG in the BS (dashed line), two DM 
models with = 500 GeV, composition W'^W~, B = 10 /rG 
and (av) = 4.3 X 10“^^ cm® s“^ in the DME region (dot-dashed 
line) and in the DMW region (three dots-dashed line). 


also note that the emission produced in the BS, in addition 
to have a weaker central emission w.r.t. the DM peaks, has 
a limited spatial extension, that does not explain all the ex¬ 
tended emission observed, especially on the north and on the 
south of this region. So, to describe the full size of the radio 
emission observed in this region, it is probably necessary to 
calculate the emission in a wider region surrounding the BS, 
or to consider the contributions coming from other sources, 
like the galaxies halos. 





Radio emission in the Bullet cluster 9 


3.3 X-ray and gamma-ray emission 

In Fig llOl we show the X-ray and gamma-ray fluxes pro¬ 
duced in the four regions considered for the models shown in 
the Fig|8] These emissions are the sum of ICS, non-thermal 
bremsstrahlung and tt'^ decay. As we can see in this figure, 
we expect a main peak of the spectral energy distribution 
(SED) around u ~ 10^® Hz coming from hadronic models 
(i.e. the WR model) and a lower one around v ~ 10^® Hz 
coming from DM annihilation models. The model is domi¬ 
nated by the emission from the MS. 

The calculated fluxes are lower than the current upper 
limits (9.1 X 10“^® erg cm“^ s“^ at iz 2.4 x 10^^ 
Hz by Fermi, and 1.1 x 10“^^ erg cm“^ s“^ at 
V - (1.2 H- 2.4) x 10^® Hz by NuSTAR). We also 

show for the sake of comparison the sensitivity curves 
for Astro-H with 100 ks of time integration (from 
http://astro-h.isas.jaxa.jp/researchers/sim/sensitivity.html), 
Fermi-LAT for 10 yrs, and CTA for 1000 hrs (from Funk & 
Hinton 2013). It is clear that the high energy emission is 
not expected to be detected with these instruments, due to 
the relatively large distance of this cluster. 

We conclude that the main source of information about 
the origin of DM and non-thermal plasmas in the Bullet clus¬ 
ter, also by considering the next generation of instruments, 
is the radio band where the very high sensitivity of instru¬ 
ments like the SKA can be used to probe the physics of this 
DM-dominated system. Note that the SKA will have also 
the possibility to measure simultaneously the magnetic field 
in the various regions of the Bullet cluster using quite pre¬ 
cise Faraday rotation measurements (see Johnston-Hollitt 
et al. 2014) and thus it will be able to fully disentangle the 
degeneracy between DM-produced secondary electrons and 
magnetic field in the various sub halos of this cluster. 
Another possibility that we discuss in next section is to use 
the SZE, by using the fact that it does not depend on the 
distance of the cluster, so it is suitable to study the Bullet 
cluster. We notice that another important property of the 
SZE is that it does not depend on the magnetic field, so it 
can allow to break the degeneracy between electrons density 
and magnetic field intensity existing in the radio continuum 
synchrotron emission. 


4 THE SUNYAEV-ZEL’DOVICH EFFECT IN 
THE BULLET CLUSTER 

The SZE is produced by the ICS of CMB photons off the 
electrons contained in ionized halos of galaxy clusters and 
other cosmic structures (Zel’dovich & Sunyaev 1969, Sun- 
yaev & Zel’dovich 1972). The SZE can be calculated for 
both thermal or non-thermal halos when the full relativis¬ 
tic formalism is used (see, e.g., Colafrancesco et al. 2003 for 
details). 

The SZE from the Bullet cluster has been already studied 
in Colafrancesco et al. (2011a) by using its spectral proper¬ 
ties over a wide frequency range from 18 GHz to 850 GHz. 
Here we perform a more detailed analysis regarding both 
the spectral and spatial properties, and using the results 
obtained in previous sections to constrain the properties of 
the non-relativistic electrons in the Bullet cluster. 
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Figure 10. X-rays and gamma rays emission in the different 
regions of the Bullet cluster with: a WR model with Sp = 2.7 
and B = 7.5 in the MS (three dots-dashed line), a WR model 
with Sp = 2.7 and B = 27 fiG in the BS (solid line), two DM 
models with = 500 GeV, composition W'^W~, B = 10 pG 
and (av) = 4.3 X 10“^^ cm® s“® in the DME region (dashed 
line) and in the DMW region (dot-dashed line). The long dashed 
line is the sum of all the components. We show also the sensi¬ 
tivity curves for Astro-H with 100 ks of time integration (from 
http: / / astro-h.isas. j axa.jp/researchers/sim/sensitivity, html), 
Fermi-LAT for 10 yrs, and CTA for 1000 hrs (from Funk & 
Hinton 2013). 


4.1 Spectral analysis 

The SZE in the Bullet cluster has been measured over a 
wide range of frequencies, with ACBAR at 150 and 275 GHz 
(Gomez et al. 2004), with the SEST telescope at 150 GHz 
(Andreani et al. 1999), with APEX at 150 GHz (Halverson 
et al. 2009), with the SPT at 150 GHz (Plagge et al. 2010), 
with ATCA at 18 GHz (Malu et al. 2010), and with Herschel- 
SPIRE at 600, 850 and 1200 GHz (Zemcov et al. 2010). 
The combination of the measurements from low to high fre¬ 
quencies allowed us to obtain information on the presence 
of multiple components in the overall SZE signal, favoring 
the presence of a second component of non-thermal origin, 
although a second thermal component with very high tem¬ 
perature cannot be excluded (Colafrancesco et al. 2011a). 

Since the SZE signal seems to be located very close to 
the MS region (see next subsection for a detailed discussion), 
we perform here a more refined spectral analysis, by taking 
into account the properties of the thermal and non-thermal 
components in the MS as discussed in the previous Sections. 

In Colafrancesco et al. (2011a), the two best possibili¬ 
ties to explain the SZE data at all the frequencies were: i) a 
thermal component with kT = 13.9 keV and optical depth 
r = 1.1 X 10“^ and a non-thermal component with a power 
law electrons spectrum with spectral index s = 2.7, normal¬ 
ized minimum momentum pi = 1 and t = 2.3 x 10““^; ii) a 
thermal component with kT = 13.9 keV and r = 3.5 x 10“® 
and a second thermal component with kT = 25 keV and 
r = 5.5 X 10~®. The first possibility was slightly favoured 
w.r.t. the second one. 

Now we can perform a more accurate analysis by setting the 
parameters of the first thermal component to the values as¬ 
sumed for the MS, corresponding to a temperature of 14.2 
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Figure 11. SZE in the Bullet cluster, with the thermal effect in 
the MS with kT = 14.2 keV and r = 1.1 X 10“^ (solid line), a 
non-thermal effect with s = 3.7, pi = 1 and t = 3 X 10“^ (dashed 
line), and the sum of the two components (thick line). Data at 
150 and 275 GHz are taken with ACBAR (Gomez et al. 2004), 
and data at 600 and 857 GHz with Herschel-SPIRE (Zemcov et 
al. 2010). 

keV (Wik et al. 2014) and an optical depth of r = 1.1 x 10~^ 
(calculated with the parameters of Ota & Mitsuda 2004), 
and considering the two cases of a second component of non- 
thermal or thermal origin. For the non-thermal case, we set 
s = 3.7, from the value of the radio spectrum in the MS (re¬ 
gion ^ 3 in Liang et al. 2000), pi = 1 (from Colafrancesco 
et al. 2011 a), and leaving r as a free parameter; for the sec¬ 
ond thermal component case we set kT = 25 keV (from 
Colafrancesco et al. 2011a), and leave r as a free parameter. 
The results are reported in the Fig llll for the non-thermal 
case, and in the Fig |12l for the two thermal components case. 
As we can see, the case with a non-thermal component pro¬ 
vides a better fit to the data than the case with two ther¬ 
mal component, especially considering the point at 857 GHz, 
hence confirming the results obtained by Colafrancesco et al. 
( 2011 a). 

It is interesting to compare the non-thermal compo¬ 
nent resulting from the analysis of the SZE spectrum with 
the electrons coming from the WR model in the MS. From 
the SZE we derived a very good agreement with the data 
by using pi = 1 , so a very low value of the minimum mo¬ 
mentum of the electron distribution is required. In the sec¬ 
ondary electron models (of both hadronic and DM origin) 
the spectrum of the electrons flattens for p ^ 10 ^, because 
of the Coulombian losses effect and of the intrinsic flatten¬ 
ing of the source term (see, e.g., Marchegiani et al. 2007 
for the hadronic case, and Colafrancesco et al. 2006 for the 
DM case), and this fact strongly suppresses the presence of 
electrons with small momentum p. So, these results point 
towards a different origin of the electrons, i.e. “primary” 
electrons, recently accelerated by some mechanism (possi¬ 
bly connected to the recent merging of the two subclusters). 
If we assume that these electrons have the same spatial dis¬ 
tribution of the thermal gas, and if we write their spectrum 
as Ns{'y) = with s = 3.7, from the value of the optical 

depth derived from the SZE analysis we obtain a normaliza¬ 
tion factor of ko = 1.0 x 10“® cm~®. In Fig ll3l we compare 
this spectrum to the spectrum of the secondary electrons de- 


Figure 12. SZE in the Bullet cluster, with the thermal effect 
in the MS with kT = 14.2 keV and t = 1.1 X 10“^ (solid line), 
a second thermal effect with kT = 25 keV and t = 1.5 X 10“^ 
(dashed line), and the sum of the two components (thick line). 
Data at 150 and 275 GHz are taken with ACBAR (Gomez et 
al. 2004), and data at 600 and 857 GHz with Herschel-SPIRE 
(Zemcov et al. 2010). 



Figure 13. Secondary electrons spectrum derived in the WR 
model in the MS (solid line) compared with the primary electrons 
spectrum derived from the SZE spectral study (dashed line). 


rived from the WR model in the MS. It is interesting to note 
that the two spectra overlap almost exactly for 7 ^ 10 ®, i.e. 
in the energy region of the electrons that produce the radio 
emission. This means that all the considerations we made for 
the WR model in the MS hold also for this model (actually, 
a slightly higher value of the magnetic field, B ~ 8 pG, is re¬ 
quired to fit the radio halo), with the difference that in this 
model the contribution of the tt® decay to the gamma-ray 
emission is not present, and therefore a much lower gamma- 
ray emission (produced only by non-thermal bremsstrahlung 
and ICS processes) is expected. We can then conclude that 
radio measurements alone cannot discriminate between the 
primary and secondary models, while the SZE observations 
indicate that a primary origin is favoured. 
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Figure 14. SZE in the Bullet Cluster produced by different com¬ 
ponents: a thermal effect in MS (solid line), a thermal effect in BS 
(dashed line) and a DM model in the DME region (dot-dashed 
line), that is the dominant component among the non-thermal 
effects in the other regions. The experimental data are not shown i 
here because no one of these models alone is able to reproduce ^ 
them (see Figs llll and [T2J. 


4.2 Spatial analysis 

There is some difficulty in the comparison between different 
maps of the SZE obtained so far for the Bullet Cluster: in 
fact, while Halverson et al. (2009) at 150 GHz and Zecomv 
et al. (2010) at 600 GHz found that the peak of the SZE 
emission is located near the MS, a bit elongated towards 
the direction of the BS, Malu et al. (2010) at 18 GHz and 
Plagge et al. (2010) at 150 GHz found that the peak of the 
SZE is shifted from the MS in the direction of the DME. 

To compare our results with these maps, we calculate 
the SZE produced in the different regions of the cluster for 
the different models used in the previous sections. There¬ 
fore, we calculate the thermal effects in the MS and BS, the 
effects for the WR models in MS and BS, and the effects 
for the DM models in DME and DMW regions (using the 
DME model normalized to the A(L) source), by using the 
full relativistic formalism (Golafrancesco et al. 2003) for all 
the models, and the formalism for the cool-core clusters for 
the BS (Golafrancesco & Marchegiani 2010). First, we com¬ 
pare the different curves so far obtained by putting them 
in the same plot (see Fig |14ll . We find that the thermal ef¬ 
fects are dominant, and the effect in the MS is dominant 
w.r.t. the effect in the BS, while the non-thermal effects by 
secondary electrons (both of hadronic and DM origin) are 
negligible, hence confirming that a primary electron pop¬ 
ulation is required in order to reproduce the experimental 
data. 

To compare the relative contributions of the secondary 
electron models in the different regions, we show in Fig |15l 
the SZE produced in the DM and WR models. The DM 
contribution results to be a factor ~ 100 smaller than the 
thermal one, while the WR contribution is smaller by a fac¬ 
tor of ~ 10®. 

In the light of these results, we conclude that the only 
important contributions come from the thermal components 
in the MS and in the BS, with the addition of the primary 
electron non-thermal effect found in the spectral analysis 
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Figure 15. Upper panel: non-thermal SZE in the Bullet Cluster 
produced in DM models in DME region (solid line) and in the 
DMW region (long dashed line), compared with WR models in 
MS and BS (dashed and dot-dashed lines), coincident with zero 
in this scale. Lower panel: non-thermal SZE in the Bullet Cluster 
produced in WR models in MS region (dashed line) and in the 
BS region (solid line). 


of the MS. This result is in agreement with the maps of 
Halverson et al. (2009) and Zemcov et al. (2010), where the 
SZE is located in the MS, with a small elongation towards 
the BS, and disagrees with the maps of Malu et al. (2010) 
and Plagge et al. (2010), where the SZE peak is shifted to¬ 
wards the DME region. We notice that a similar elongation 
seems to be present also in the radio maps (Shimwell et al. 
2014), and that in this region there is a high temperature 
region with temperature higher than 20 keV (Markevitch et 
al. 2002). This last point can explain why the SZE appears 
to be shifted towards this region, being the SZE propor¬ 
tional to the Gompton parameter y oc J dlnkT, while in 
this region the X-ray emission is low, so the density of the 
gas should be small (and it does not emerge also at larger 
X-ray energies, see Wik et al. 2014). Since in this region it 
is present the radio source A(L) that we already studied for 
its radio emission, it is also possible that this galaxy yields a 
contribution to the total SZE (see, e.g., Golafrancesco 2008, 
and Golafrancesco et al. 2013 for a study of the SZE in the 
radio galaxies lobes). 

In Fig |16l we show some simulated maps of the SZE at 
several frequencies produced with the different models found 
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in the previous sections, including for the MS the contribn- 
tion from the non-thermal SZE from primary electrons. 

At 18 GHz (the frequency of the map of Main et al. 
2010) the SZE is very low, because this is a small frequency 
w.r.t. the frequency of the minimum; the SZE is negative 
and is dominated by the thermal effect in the MS, with 
a small elongation towards the BS and two point sources 
corresponding to the DM peaks. 

At 150 GHz all the effects are negative (this frequency 
is near to the minimum of the SZE), and it is interesting to 
note that the combined effect of the SZE from MS and BS 
produces a secondary SZE peak in the middle of these two 
regions. This can be the origin of the elongated shape found 
in the map of Halverson et al. (2009) at the same frequency. 
The DM peaks appear to be very narrow, and they can be 
considered as point-like sources for most SZE experiments. 

At 217 GHz, the frequency of the crossover point of 
the thermal effect in non-relativistic limit, the SZ effects are 
still negative and for the MS the dominant component is 
the non-thermal one for primary electrons, while there is 
still a residual thermal effect from the BS. We show also the 
case for 224.5 GHz, very close to the crossover frequency 
for the thermal effects if relativistic effects are considered, 
where the BS disappears almost entirely, and only the non- 
thermal component in the MS is present. This can be a good 
frequency where it is possible to test this model by detecting 
the non-thermal SZE (if we can distinguish it from a possible 
kinematic effect). 

At 450 GHz the map is dominated by the thermal effects 
in the MS and BS, with the peak in the middle present also 
in this case. We note that, while the overall effect is positive, 
the DM peaks are negative, so they appear like two small 
holes in the map. 

At 600 GHz (the frequency of the map of Zemcov et 
al. 2010) the SZE is positive, and it is dominated by the 
thermal peaks. 

At 850 GHz the non-thermal effect in the MS is dom¬ 
inating w.r.t. the thermal effect, and the thermal effects in 
MS and BS appear as small corrections; the DM components 
are still negative and negligible w.r.t. the other effects. 

We can compare these results with the expected perfor¬ 
mance of the Millimetron mission. Millimetron is planned 
to have, at the wavelength of 300 fim and when operating 
in the single dish and medium spectral resolution mode, an 
angular resolution of « 6 arcsec and a sensitivity of « 4 pjy 
for 3600 seconds of integration (Kardashev et al. 2015). In 
the SZE maps of Fig |16l the lowest level contour is 1 mjy 
arcmin“^; the flux integrated inside a circle with diameter 
of 6 arcsec for the lowest level contour of our maps is ~ 8 
/rjy, larger than the previous sensitivity limit. For this rea¬ 
son, we expect that Millimetron will be able to check our 
predictions by detecting the flux down to the lowest contour 
level we show in these maps, allowing hence to better con¬ 
strain the properties of the different subhalos of the Bullet 
cluster. 


5 DISCUSSION AND CONCLUSIONS 

In this paper, we studied the complex structure of the radio 
halo in the Bullet Cluster by considering the properties of 
each of the baryonic and the DM subhalos observed with 


X-rays and gravitational lensing measurements. The radio 
emission from the MS is well explained by assuming a mag¬ 
netic field of ~ 7.5 /rG, and a population of relativistic elec¬ 
trons centered on the MS. The radio observations do not 
allow to discriminate between a primary and a secondary 
origin for these electrons, and the forthcoming instruments 
in the X-ray and gamma-ray bands like Astro-H and CTA do 
not seem to be able to solve this problem, because the high 
distance of this cluster causes the emissions in these spec¬ 
tral bands to be smaller than the sensitivity limits of these 
instruments. Instead, SZE observations at high frequencies 
indicate that a primary origin is favoured. 

The study of the other subhalos from radio data is more 
complicated, mainly because the available data do not sep¬ 
arate the emissions coming from each of them, but consider 
the emission from a big region including all the subhalos and 
several point and extended sources. The radio spectrum of 
this region indicates that the presence of different emissions 
with different spectra is very probable. The combination of 
the emissions coming from the BS and the DME region can 
produce a spectrum similar to the observed one, but there 
are still various open issues: 

i) the curvature of the radio spectrum at v > 5 GHz is dif¬ 
ficult to be obtained for all the DM models we considered; 
this problem could be less important if we think that the 
procedure of sources removal in this region can be delicate; 
a) the required annihilation cross section is very high 
{(crv) ~ 8.1 X cm® s“®), a factor of order of ~ 10® 

larger than the present upper limits, and probably it cannot 
be reproduced by a boosting factor due to DM substructures 
(see, e.g., Colafrancesco et al. 2015); anyway, if a higher (but 
reasonable) magnetic field of 5 pG is assumed, a similar spec¬ 
trum can be obtained with a cross section larger by a factor 
of 10® w.r.t. the present upper limits, and this difference can 
be recovered with a DM substructure boosting factor; 

Hi) for these DM properties, an emission coming from the 
DMW region should also be observed, and this does not 
happen. 

To better constrain the DM properties it would be 
important to perform better observations of each of the 
radio subhalos separately; this aim, both considering the 
sensitivity and the angular resolution properties, can be 
reached by forthcoming instruments like SKA, or its pre¬ 
cursor MeerKAT. 

Another possibility is that the radio emission coming 
from the DME region is really visible as a point source (due 
to the central peak of the DM distribution) that has been 
confused with the emission of the A(L) radio source, while 
its extended halo is included in the extended radio emission 
observed in this region. This hypothesis requires a smaller 
value of the cross section, and therefore also a smaller radio 
emission from the DMW region. To test this hypothesis, very 
detailed radio measurements, especially near the DME and 
DMW regions, are required. 

A different possible explanation of the radio emission 
close to the DME region can be obtained by considering 
that near the DM peak there is also a big concentration of 
galaxies, so it is possible that the extended radio emission 
observed in this region is produced by the electrons emitted 
by the galaxies and diffusing in the IGM (Colafrancesco et 
ah, in preparation). This is probably true at least for the 
emission coming from the regions in the north and the south 
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Figure 16. First panel: radio map at 1.5 GHz simulated with the surface brightness profiles plotted in the Fig|9] Other panels: SZE 
maps at different frequencies simulated with the same models, and adding a non-thermal component in the MS produced by primary 
electrons as found in Sect.4.1. In all the panels, on the axes are the values of the coordinates differences (in arcmin) w.r.t. the origin, 
fixed in the point 06:58:50 -55:59:00 (J2000). In the first panel, contour levels correspond to: (1 X 10“®, 5 X 10“®, 1 X 10“®, 5 X 10~®, 
1 X lO-"', 5 X 10"'^, 1 X 10-3, 5 X 10-3, 1 X 10-2, 5 x IQ-^, 1 x IQ-l, 5 x IQ-l, 1, 5, 10, 50 and 100) Jy arcmin-^. In the other panels, 
contour levels correspond (in absolute value) to: (1, 5, 10, 15, 20, 25, 30, 35) mjy arcmin—2. 


of the BS, where the small extension of this subcluster does 
not predict a strong emission that instead is observed. We 
will test this hypothesis in a future paper. 

Finally, the SZE appears to be the best tool to constrain 
the properties of the relativistic electrons in addition to the 
radio measurements. While the SZE produced by electrons 
produced by DM annihilation is probably too small to be 
observed, the study of the non-thermal SZE is anyway im¬ 
portant in order to determine the properties of the electrons 
of baryonic origin, and to remove their contributions from 


the overall emission, favoring the study of the properties of 
the DM. Future instruments like Millimetron will be crucial 
to this aim. 
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